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The MILC Collaboration has been producing ensembles of lattice configurations with three dynamical flavors 
for the past few years. There are now results for three lattice spacings for a variety of light and strange quark 
masses, ranging down to mi = 0.1m s , where m s is the dynamical strange quark mass and m; is the common mass 
of the u and d quarks. Recently, the Fermilab, HPQCD, MILC and UKQCD collaborations have presented a 
summary of results obtained using these lattices. Compared with quenched results, these new calculations show 
great improvement in agreement with experiment. This talk addresses the technical improvements that make 
these calculations possible and provides additional details of calculations not presented in the initial summary. 
We demonstrate that a wide range of hadronic observables can now be calculated to 2-3% accuracy. 



1. INTRODUCTION 

The past several years have seen a marked im- 
provement in our ability to do realistic calcu- 
lations in Lattice QCD (LQCD). Improved ac- 
tions have better scaling properties thus allow- 
ing us to more quickly approach the continuum 
limit than with the simplest actions. For stag- 
gered quarks, these improved actions also allow 
us to perform simulations (on current comput- 
ers) with small enough up and down quark masses 
that we can use chiral perturbation theory to ex- 
trapolate to the physical quark masses. Further, 
the use of improved actions has greatly reduced 
the "flavor" symmetry violations, that we now 
call "taste" symmetry violations. Without this 
improvement, many of the non-Goldstone pion 
states are as heavy as the lightest kaon state. 
However, with the improvement in taste symme- 
try, the splitting between pions and kaons is more 
clear, and one expects to be better able to discern 
the effects of a dynamical strange quark. These 
improvements, and the general increase in avail- 
able computer power have made it attractive to 
embark on a large scale calculation of QCD with 
three dynamical quarks. 

Generation of the gauge ensembles has been 
going on for four years. MILC has been shar- 
ing dynamical configurations through the Gauge 
Connection and informally for some time. The 



Fermilab, HPQCD, MILC and UKQCD collabo- 
rations have all been using the MILC configura- 
tions to calculate a variety of quantities with both 
light (u, d or s) and heavy (c or b) quarks. Taken 
together, we find that these calculations result in 
a compelling picture. LQCD with three flavors 
agrees with experiment with errors of 2-3% for 
a wide variety of quantities, whereas quenched 
QCD has errors as large as 15-20%. 

In the rest of this article, I shall briefly discuss 
some of the issues that arise when dealing with 
staggered quarks. I will then give an introduc- 
tion to the parameters set in lattice calculations, 
and some details of the MILC ensembles. After 
an explanation of high precision observables, we 
go on to discuss a variety of calculations such as 
the light meson decay constants, the heavy quark 
spectrum, the light quark spectrum and topolog- 
ical susceptibility. We conclude with a discussion 
of future prospects. 

2. ALGORITHM ISSUES 

Staggered quarks are fast compared to other 
discretizations, but one must take a fourth root to 
deal with the so-called fermion doubling problem. 
This leads to a potential loss of locality, and we 
don't know how to construct a transfer matrix. 
Nevertheless, there are positive things one can 
say. 
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If one considers a formal perturbative expan- 
sion of the theory, the fractional root causes no 
problems in the expansion. 

It is known that the CP violation that should 
occur when m u + m,4_ < does not happen with 
staggered quarks, but this deficiency is not obvi- 
ously relevant to the real world. 

The taste singlet axial current is only partially 
conserved so the chiral anomaly is preserved. 

We now understand how to modify chiral per- 
turbation theory in the presence of taste changing 
interactions jll2j . 

An important part of recent progress in LQCD 
is the use of improved actions. The MILC collab- 
oration uses the "Asqtad" action [3I4I5I6I7] . 

The gauge action includes the plaquette, the 
1x2 rectangle and a bent parallelogram 6-link 
term. The quark action includes paths up to 
seven links long as shown below and the 3-link 
Naik term. 

Each diagram in Fig.^represents a term in the 
quark action of the form x(x)V(x,x + fi)x(x + fi), 
where V is the product of links along the path. 

The "fat-link" contributions help to suppress 
taste symmetry breaking. The final 5-link path 
was introduced by Lepage jS] to correct the small 
momentum form factor. The result is an action 
that has leading errors of order a 2 a and a 4 . 

The decision to base our calculations on the 
Asqtad action is based not on a fondness for frac- 
tional powers, but on the expectation that we can 
get to the correct chiral physics more easily this 
way. The chiral extrapolation for the u and d 
quarks has been a challenge for LQCD through- 
out its history. 

To generate ensembles of gauge configurations, 
we must select certain physical parameters: the 
lattice spacing (a) or gauge coupling (/?) , the grid 




Figure 1. Some of the products of links used in 
the Asqtad action. The Naik term is not shown. 



Table 1 

Masses, coupling and number of configurations in 
MILC ensembles 



a = 0.12 fm; 20 3 x 64 (coarse) 


am Ut d / am s 


10/.g 2 


# conf. 


0.40 /0.40 


7.35 


332 


0.20 /0.20 


7.15 


341 


0.10 /0.10 


6.96 


340 


0.05 /0.05 


6.85 


425 


0.04 /0.05 


6.83 


351 


0.03 /0.05 


6.81 


564 


0.02 /0.05 


6.79 


485 


0.01 /0.05 


6.76 


608 


0.007/0.05 


6.76 


447 


0.005/0.05 


6.76 


137 


a = 0.09 fm; 28 3 x 96 (fine) 


am Utd 1 am s 


lO/.g 2 


# conf. 


0.031 / 0.031 


7.18 


336 


0.0124 / 0.031 


7.11 


531 


0.0062 / 0.031 


7.09 


583 



size (-/V 3 x N t ), and the quark masses (m„ i£ j, fn s ). 

To control systematic errors, we must take the 
continuum limit, take the infinite volume limit 
and extrapolate to light quark mass for the u and 
d quarks. However, we can work at the physical 
s quark mass. 

In the Table ^ we give some of the details of 
the MILC ensembles. We are able to go closer 
to the chiral limit than we did with unimproved 
staggered quarks in the mid-90s. This can be seen 
in Fig. |21 where the pion and nucleon masses are 
plotted in lattice units. The straight line corre- 
sponds to the physical pion to nucleon mass ratio. 
Approaching this line corresponds to the chiral 
limit. As the plot is a log-log plot, approaching 
the lower left corner (and going beyond) corre- 
sponds to the continuum limit. Note the three 
octagons to the left of the lowest left-most dia- 
mond. At this lattice spacing, with 2+1 flavors 
we have gone closer to the chiral limit. One also 
sees fancy plus symbols at smaller lattice spacing 
and comparable or smaller m^/mjv. 
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□,x quenched 
n,o,w three flavor 
o unimp. two flavor 




n F — results 



n r =3 results (2003) 



Figure 2. Nucleon and pion masses for MILC 
unimproved two-flavor staggered runs and im- 
proved runs with Nf = or 3. 



3. RESULTS FOR HIGH PRECISION 
OBSERVABLES 

Certain quantities are intrinsically easier to cal- 
culate on the lattice than others. For example, a 
ground state mass is much easier to calculate than 
an excited state mass in the same channel. Also, a 
stable particle mass is much easier to calculate ac- 
curately than an unstable particle mass. The Fer- 
milab, HPQCD, MILC and UKQCD collabora- 
tions have been working together to confront a va- 
riety of experimental results, and to make predic- 
tions that can be tested. So far, we have mainly 
been concentrating on high precision observables. 
If we fail to reproduce well known experimental 
results, after having controlled the chiral, contin- 
uum and finite volume extrapolations, we are left 
with little wiggle room. That is, we can no longer 
blame the quenched approximation. We would 
either conclude that staggered quarks are flawed, 
or that there is a problem with (L)QCD. On the 
other hand, if we are able to demonstrate good 
agreement with a variety of measured quantities, 
it would give us confidence that we are also able 
to make predictions, many of which are crucial 
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Figure 3. Ratio of lattice result to experimental 
value for nine quantities. Quenched (dynamical) 
results on the left (right). 



for accurate determination of standard model pa- 
rameters such as quark masses and CKM mixing 
matrix elements. 

Our results depend on four input masses and 
the gauge coupling (equivalently the lattice spac- 
ing). The lattice spacing is determined from the 
T' — T mass difference, a quantity with little de- 
pendence on the b quark mass [Sj. The common 
mass of the u and d quarks can be determined 
from TO^, while m s is determined from 2m^-- to J, 
or a simultaneous fit to m n and ttik determines 
m Ut d and m s . The charm quark mass m c is de- 
termined from m p a , and is determined from 
mx- 

The results that follow are mainly the work of 
the following people and collaborations [§]: 

C. T. H. Davies, E. Follana, A. Gray, G. P. Lep- 
age, Q. Mason, M. Nobes, J. Shigemitsu, 
H. D. Trottier, M. Wingate: HPQCD and 
UKQCD Collaborations; 

C. Aubin, C. Bernard, T. Burch, C. DeTar, 
S. G., E. B. Gregory, U. M. Heller, J. E. Het- 
rick, J. Osborn, R. Sugar, D. Toussaint: MILC 
Collaboration; 

M. Di Pierro, A. El-Khadra, A. S. Kron- 
feld, P. B. Mackenzie, D. Menscher, J. Simone: 
HPQCD and Fcrmilab Collaborations. 

The work on the topological susceptibility was 
done by MILC including T. DeGrand, A. Hasen- 
fratz and A. Hart [TUj. 
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4. 7T, K MASSES AND DECAY CON- 
STANTS 

In our latest calculations, we have a greatly im- 
proved understanding of chiral extrapolation due 
to five factors: 1) inclusion of taste symmetry 
breaking, 2) improved action, 3) more and lighter 
quark masses, 4) partial quenching and 5) simul- 
taneous fit to mass and decay constant results. 

In MILC's prior work with two flavors of unim- 
proved KS quarks, we had five dynamical quark 
masses at several fixed values of j3. With no con- 
sideration for taste symmetry breaking, we tried 
seven combinations of powers and logs to try to fit 
the pion masses, but our best fits had only a com- 
bined confidence level of 0.1. In that work, the 
valence light quark mass was the same as the dy- 
namical mass, so there were only five data points 
for each (3 value. (In the current work with par- 
tial quenching, we have over 100 points for each 
lattice spacing.) 

In Fig. 2] we plot some of our data for the five 
ensembles with the lightest m u> d on the coarse 
lattice JT]. Quark propagators for nine differ- 
ent valence quark masses were calculated. The 
octagons correspond to equal mass quark and an- 
tiquark. The crosses have the antiquark fixed at 
the dynamical strange quark mass. Many addi- 
tional points with different combinations of the 
valence masses are not plotted. The bursts in- 
clude results from ensembles with heavier quarks 
than those listed in the legend. 

The leading order formula from chiral per- 
turbation theory for pion mass and decay con- 
stant depends on four combinations of the Gasser- 
Leutwyler constants. 
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Figure 4. Some of trie" resufls 1 for 7'ito 2 /(toi + 
m?) on some of the coarse ensembles. There are 
too many partially quenched results to plot all 
available data. 



However, with staggered quarks there are taste 
symmetry violations that modify these formu- 
lae to take into account the fact that all the 
pion states are not degenerate with the Goldstone 
pion. These formulae are quite complicated and 
may be found in Ref. |12|. They involve new cou- 
plings 8y and S' A from the taste symmetry break- 
ing part of the effective Lagrangian. The formulae 
also involve the chiral log function, i, given by 

,2 



TYl 

In —j- + 5\{mL) 
A 

x 



(3) 



(2) 



with A x the chiral scale, 8\ the known finite vol- 
ume correction and L the spatial size. 

In Figs. and we show results for 
rim% r /(m x + m y ) and f n ri with m x and m y the 
two valence quark masses, and n a scale set 
from the heavy-quark potential. Only some of 
the data points are plotted. The curves corre- 
spond to partially quenched results on ensembles 
with the light quark mass indicated in the legend. 
Coarse (fine) lattice results are plotted with dia- 
monds (squares). The confidence level of the fit 
is 0.72, ignoring the contribution of the Bayesian 
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CL = 0.72; 0.73 
o a=0.13 fm; □ a=0.09 fm 

m L =0.6md 
0.4m ri 




physical m u + m d 



0.00 0.05 0.10 0.15 

(m x +m )r, 

Figure 5. Selected data for the pion mass and fit 

curves from simultaneous fit to mass and decay 
constant results. 
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Figure 6. Selected data for pion decay constant 

and fit curves from simultaneous fit to mass and 
decay constant results. 



constraints and 0.73 including them. The points 
are plotted after taking into account finite vol- 
ume corrections that are sometimes as large as 
1%. In Fig. [5] a vertical bar shows the value of 
m u +rrid, determined from the physical pion mass. 
The extrapolated value for f n and the experimen- 
tal result are shown in the lower left of Fig. HjJ 
Note the close agreement. The systematic error 
includes variations in which points are included 
in the fit, variations in assumptions about how 
rapidly the taste-violating terms and non-taste 
violating terms are changing with lattice spacing, 
and, most important, the scale error of 2.2%. In 
Fig. [7] a similar plot is shown for fx ■ The differ- 
ence here is that more heavy valence mass points 
are included in the fx fits. One can see that 
the Bayesian constraints contribute more to \ 2 
here. Agreement with experiment is also accept- 
able. The dynamical strange quark mass used in 
these ensembles is 10-20% higher than the result 
obtained for the valence m s . Additional coarse 
ensembles are being generated with am s — 0.03. 

Figure |H1 shows the continuum limit of f n and 
fx- For these decay constants, the leading cor- 
rection is linear in a s a 2 . Values for the coarse 



and fine results are also shown with taste sym- 
metry violating terms set to zero. When this is 
done, there is very little difference between the 
two lattice spacings, supporting the contention 
that taste violations are the major source of dif- 
ference. 

4.1. Gasser-Leutwyler Parameters 

We now have preliminary results for several 
combinations of Gasser-Leutwyler parameters us- 
ing m, as the chiral scale. We find 

2L 6 -U = 0.5(2)(tJ) x 10~ 3 (4) 

2L 8 -L 5 = -0.1(1)(±|) xHT 3 (5) 

U = 0.3(3)(±£) x 10~ 3 (6) 

L 5 = 1-9(3)(±|) x 10- 3 , (7) 

where the first error is statistical and the second is 
systematic. The systematic errors are dominated 
by differences over various acceptable chiral fits. 

In the continuum, it is expected that L 5 = 
2.3(2) x 10~ 3 and L 4 « L 6 rj [T3]. 

The result for 2L$ - L 5 of -0.1(1) (1^) x 10~ 3 
is well outside the range that allows for m u = 
|14| which is (approximately) 

-3.4 x 10~ 3 < 2L$ -L 5 < -1.8 x 10~ 3 . (8) 



6 



0.22 



0.20 



s-, 



0.18 



0.16 




O m 
□ m 



0.6m a 
0.4mS 
0.2m 
L -0.14rfig' 

0.2m! 



m„=0.8ms; m 



extrap; ISyStematic err 
experiment 



0.10 



0.15 



0.20 



(m +m )r. 

Selected data for kaon dc 



Figure 7. Selected data for ka&n decay constant 
and fit curves from simultaneous fit to mass and 
decay constant results. 
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kaon decay constants. 



Thus, if MILC's preliminary result holds up, it 
would rule out the m u = possibility. This re- 
sult seems likely to be robust under changes in 
treatments of higher order terms. 

5. HEAVY QUARK SPECTRUM 

The most stringent tests of LGT come from 
hadrons that require no chiral extrapolation for 
the valence quarks, are stable to strong decay, 
and far from threshold. For example, J/^P, T, 
D s and B a meet these conditions. The HPQCD 
and UKQCD collaborations have used a non- 
relativistic QCD (NRQCD) ^H] action accurate 
to order v 4 to study bottomonium on MILC con- 
figurations. The Fermilab collaboration has used 
clover quarks with the Fermilab interpretation 
to study charmonium states. 

The spin averaged bottomonium splittings are 
shown in Fig. [5| for the IP, 2P and 3S differences 
from the IS level [T7]. (Recall, the 2S-1S splitting 
is used to set the scale.) Three coarse dynamical 
ensembles have been used. The experimental re- 
sults (without error bars) are plotted near the 
left edge of the diagram. On the right hand side 
of the diagram, the quenched values are plotted. 



The corresponding values in Fig. Income from the 
quenched results and lightest dynamical mass re- 
sults shown here. 



Several fine and hyperfine splittings are shown 
in Fig. 1101 which contrasts results in the quenched 
approximation with results on one dynamical en- 
semble with m u ^d = 0.2m s . There are significant 
differences for the hyperfine splittings of T and 
T'. However, the corresponding rj states have not 
been observed. The fine structure in the x states 
is improved on the dynamical ensemble. 



Another quantity in the ratio plot is 2m b — mr- 
Using NRQCD, this quantity has been calculated 
on three of the ensembles (m U: d = 0.01, 0.02 and 
0.03). The splitting is shown in Fig. ^2f° r var- 
ious valence quark masses. The bursts are the 
experimental values. For the ratio plot (Fig.EJ, a 
correction has been applied to account for the fact 
that the dynamical strange quark mass is about 
20% too large. 
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NRQCD. 
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Figure 10. Several hne and hyperhne splittings 
in the T system using NRQCD. Experimental 
results indicated with dashed lines; open circles 
denote quenched approximation, filled circles are 
from the coarse ensemble with am u ^d = 0.01, 
to, = 0.05. 



The Fermilab collaboration has looked at the 
spin averaged IP-IS splitting for charm quarks 
using clover/Fermilab quarks Four dynami- 
cal ensembles with am u 4 = 0.007, 0.01, 0.02 and 
0.03 have been analyzed. Results are shown in 
Fig. El O n the same ensembles, the hyperfinc 
splitting between ip and -q c has been calculated. 
It is about 80% of the experimental value. This 
splitting is sensitive to the coefficient of a^F^. 
This coefficient needs to be calculated to 1-loop 
level. It is hoped that using the corrected coeffi- 
cient will improve the accuracy of this result. 

6. LIGHT QUARK SPECTRUM 

Hadrons containing up or down valence quarks 
require chiral extrapolation before they can be 
compared with experiment. We have seen how 
this can be done including taste symmetry break- 
ing terms for the pseudoscalars. For the nucleon, 
there is not yet a corresponding formula. In 
Fig- El we show rnjvri for the coarse and fine en- 
sembles as well as some ensembles with a w 0.20 
fm. The figure contains a phcnomcnological chi- 
ral perturbation theory curve, which is not a fit 
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Figure 13. Hyperfine splitting of the IS charmo- 
nium levels in lattice units vs. the light dynam- 
ical mass in lattice units from four of the coarse 
ensembles. 



to the data. Although this figure may be sug- 
gestive that the lattice results will approach the 
experimental result in the chiral and continuum 
limit, more work is needed. The nucleon is not 
yet a high precision lattice quantity. 

For the nucleon, we need to consider finite size 
effects carefully. All our previous studies were 
done with two dynamical flavors and with much 
larger taste symmetry breaking. This case may 
have smaller finite size effects. For m Ut d — 0.2m s 
we are studying 28 3 x 64 to complement the 
20 4 x 64. Preliminary results shows about a 1% 
decrease in mass on the larger volume. More 
statistics are needed. Another critical need is a 
chiral fit to the nucleon masses including taste 
violations. 

With the additional caveat that the p can decay 
and the lattice fits take no particular account of 
this possibility, I present MILC's current APE 
plot in Fig. E>1 



7. TOPOLOGY 

At Lattice 2002 there was a concern that the 
dynamical quarks were not sufficiently suppress- 
ing topological susceptibility at light quark mass 
|19|. Since then, the fine lattice runs have been 
extended [TJ31- This has resulted in smaller errors, 
and the 3-flavor point that was so high previously 
has come down. Figure lTBI shows an extrapolation 
to the continuum limit in the quenched approx- 
imation and for the three fine lattice ensembles. 
We see in Fig. El that the continuum extrapola- 
tions are in reasonable agreement with the leading 
order theory 



8. PROSPECTS 

We feel that with the Asqtad action and 
2+1 dynamical quarks, we are making excellent 
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Figure 17. Current topological susceptibility re- 
Figure 15. APE plot for Asqtad improved KS sults Qn Asqtad ensembles . 

quarks, with 0, 2 and 3 flavors. 
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progress. Most notably, the ratio plot is much 
improved compared with the quenched approxi- 
mation. This is particularly true for the 6-quark 
spectrum. Much other work is in progress that 
could not be discussed in this talk, and much ad- 
ditional work remains: 

• Control chiral limit better on fine lattice 

• Study possible finite size effects further 

• Reduce lattice spacing to demonstrate con- 
trol of continuum limit 

• Complete physics measurements on more 
ensembles, fg, form factors, etc. [2U| 

• Use a highly improved FNAL-type heavy 
quarks [21] 

• Develop new techniques for unstable parti- 
cles and excited states 

• Complete additional perturbative calcula- 
tions for renormalization and matching |22| 

This work was supported by the U.S. De- 
partment of Energy under grant FG02-91ER 
40661. I am pleased to thank these collabora- 
tors who provided graphs or commented on the 
manuscript: C. Bernard C. Davies, C. DeTar, 
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